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enhanced because of large number of KK modes in models with right handed singlet 
neutrinos in large extra dimensions. We find that the invisible decay rate of Higgs can 



X 



In light of current atmospheric neutrino oscillation data, we revisit the invisible 
decay of the standard model Higgs boson and other pseudoscalar mesons which can be 



be as large as H — > bb decay rate only for a very restricted region of parameter space. 
This parameter space is even further restricted if one demands that the dimensionless 
neutrino Yukawa coupling A is O(l). We have also studied the scenarios where singlet 
neutrino propagate in a sub-space, which lowers the string scale M* and keeps neutrino 
Yukawa coupling 0(1). We have also considered decays of other spin-0 mesons to vv 
and found the rates to be too small for measurement. 
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The concept of large extra dimensions and TeV scale quantum gravity [1] introduced by 
ADD has attracted a lot of attentions recently. In this scenario, one has 5 additional spatial 
dimensions of size R in which gravity propagates whereas the standard model particles with 
chiral fermionic content are confined to the usual 4-dimensions (4D). The effective AD 
Planck scale, M Pi ~ 2.4 x 10 18 GeV , is related to the fundamental Planck scale M* in 
(4 + 5) dimension by 

Mp t ~ R 5 M 5 + 2 (1) 

Thus, for the large extra-dimensions, it is possible to have a fundamental scale M* as low as 
a TeV [1], helping to resolve the gauge hierarchy problem of the standard model. The size of 
the extra dimensions can be as large as ~mm for 5 = 2. Another interesting aspect of this 
scenario is the generation of a small neutrino mass [2, 3]. The relatively small value of 
indicates the possibility of a seesaw mechanism, with right-handed (RH) singlet neutrinos 
also propagating in the full 4 + 5 dimensions with gravity [4, 5, 6, 7, 8, 9]. For illustration 
let us assume the case of a single extra-dimension labeled by y (x labels the usual 4D) [10]. 

A massless Dirac fermion A" which is a standard model singlet lives in 5D. The T matrices 
in 5D can be written as 

I \ * I i \ 

r" = ,r 5 = . 2 

\ff" / \ -i J 

The Dirac spinor in the Weyl basis can be written as 

»-(*). 

where ip and x are 2-component complex spinors (with mass dimension 2). The 5D kinetic 
term for is 

%ee = / d4xd y N ( rM ^ + r5 ^) N - ( 4 ) 
and the interaction action is 

A 



S int 



= Jd*x-j= l(x)H*(x)^(x,y = 0), (5) 



where, £ = (e, v) is the standard model lepton doublet and H is the Higgs doublet and A is a 
dimensionless Yukawa coupling in 5D. We have assigned and £ opposite lepton number, 
so that the lepton number is conserved in Eq.(5). 

In the effective 4D theory, A^ appears as a tower of Kaluza-Klein (KK) states: 

i> = J2^ n \x)e in «/ R , (6) 

n V-K 



where ip^ are AD states. Similarly, x has a KK tower, x^ ■ Thus, in AD, we have the 
effective action: 

%ee = / d^T, W n) ^ {n) + X (n) ^X (n) + (^ (n) X (n) + h.c.)l , (7) 
J n L \K J . 

S int = Jd'x £ \^l{x)H*{x)^\x). (8) 



where n/R is the Dirac mass for the KK states. 

Thus, we see that the AD neutrino Yukawa coupling is suppressed by the volume of the 
extra dimensions, or in other words, by the ratio M*/Mpj. When H gets a VEV ( v ~ 246 
GeV ), we get Dirac mass for standard model neutrino (v) denoted by m: 

AM* M* 4 

ra ~ — ^— — v ~ A— —10 eV. (9) 
v^M Pi TeV V ; 

However, this standard model neutrino (which is dominantly the lightest neutrino) has a 
small mixture (~ mR/n) of heavier neutrinos. With this, the lightest neutrino mass is 
modified to ( we extend the formula to 5 extra-dimensions) 

Tfl 

m u a* = (10) 



1 + 



m 2 M z pl 271-V 2 
Mjf r(5/2)(<5-2) 



Now, for a given value of M* and 5 the upper limit of m v is m™ a:E ps -j^ \J ^ ^// 2 ^ , ( for 
5 = 2, 5 - 2 is replaced by l/ln(Mpi/M*)). We then get 

m , . 

m v pa (11) 

i+( m y y ' 

Invisible decay of Higgs 

It has been shown in Ref.fll] that the standard model Higgs can decay into v L v l R with 
the strength ~ ^Jfjr^ which looks like rather a small number. However, one can a get large 
enhancement in the decay rate when one sums the KK excitations of right-handed neutrino 
states with masses below the Higgs mass (M H ). This summation is proportional to the 
volume R 5 of the extra-dimensional space, with a momentum-space factor of order (Mh) s . 
So, the rate is enhanced by the factor (Mpi?) 5 . Therefore the sum of partial widths of the 
Higgs into KK excitations of neutrinos is of the order of : 

E r(# - vrf R ) ~ ^Y?{M H R)> (12) 



where, Y v ~ ^it^- To quantify the invisible decay rate of the Higgs boson, following the 
notation of Ref.[ll], we define the ratio k(= J2iB(H — > v L is l R )/B(H — > bb)), which can be 
expressed as : 

rn^(MH\' (Mpi\ 2 {m 

From this expression it is clear that k increases with m and M H , while it decreases with 
increase in 5 and M*. Note that k depends on m rather than m v directly. 

We now evaluate the ratio k in light of the latest atmospheric neutrino oscillation data, 
which suggest Am 2 tm = (1.5 — 4.0) x 1CT 3 eV 2 [12]. In our analysis we fix the highest 
neutrino mass which dominates H — > vv rate to be m v ~ \j Aml tm and then solve for m from 
Eq. (10) in terms of M* and 5. The seesaw mechanism in Eq. (9) clearly favors hierarchical 
neutrino masses rather than degenerate mass spectrum ( which will need fine tuning of A). 
We consider two cases that respect the allowed range of Am^ m , (a) m? v = 1.5 xlO" 3 eV 2 and 
(b) ml = 4.0 x 10~ 3 eV 2 . Another important consideration is the value of A that ensues for a 
given m and M*. Clearly, the value of A, which is the higher dimensional Yukawa coupling, 
should not exceed ~ 10 — 20, otherwise we will be in a non perturbative regime. This is a 
serious constraint on allowed value of M* for a given 5. 

In Figure 1 (a) we show the variation of k with M* for three choices of extra dimensions 
5 = 3,4 and 5 keeping m 2 v = 1.5 x 10~ 3 . The solid, broken and dotted curves correspond to 
100, 150 and 200 GeV Higgs mass respectively. Left end of each curve has been truncated 
by using a constraint m < 4m u . Value of m higher than this correspond to unacceptably 
large A. From this Figure, it is clear that for 5 = 4 and 5 the invisible width is only a tiny 
fraction of H — > bb rate. For 5 = 3, the value of k is at most 1, 3 or 7 for the Higgs masses 
100, 150 and 200 GeV respectively, and less than 1% for most of the parameter space. 

Similar variation is shown in Figure 1 (b), however, for ml ~ 4.0 x 10~ 3 . The situation 
is better in this case, since, higher value of m v lead to larger value of m in the ratio k. In 
this case also we have truncated each curve on the left by imposing the same condition as 
in Figure 1(a). 

Before we proceed, we would like mention two points here. Firstly, beyond M H = 150 
GeV, the dominant mode of Higgs decay is if — ^ WW*, thus the H — > bb branching ratio 
is very small, and any H — > uu branching ratio will also become small. However, for the 
purpose of comparison we have presented the case of 200 GeV Higgs mass. One should also 
note that for the value of M* at which k is 0(1) and larger, corresponding Yukawa coupling 
A is around 80 for 8 = 3. Such a large value of A may not be accepted from the perturbative 
point of view. 



For a given value of m v and M*, the value of A can be obtained from Eq. (9) and Eq. 
(10): 



A = {mjMj 10 16 

Mj~Mj" r(<5/2)(<5-2) 



It is easy to show that if we limit A < 10 (to be in perturbative regime) , for m 2 = 1.5 x 
10~ 3 eV 2 we are restricted to M„ > 40 TeV and for m 2 = 4.0 x 10~ 3 eV 2 , M* > 64 TeV. The 
dependence on 5 is very week on these bounds. Thus for such high values of M* the ratio k 
is highly suppressed. 

For completeness we discuss the case of degenerate neutrinos, separated by Am 2 consis- 
tent with the atmospheric and solar neutrino data. Recently WMAP [13] provided important 
new information on cosmic microwave background anisotropies. After combining the data 
from 2dF Galaxy Redshift Survey, CBI, and ACBAR[14], WMAP places stringent limits on 
the contributions of neutrinos to the energy density of the universe [15] 



which implies 



Q u h 2 = < 0.0076 (95%C.L.), (15) 

93.5 ev 



<0.71eV (16) 



for a single active neutrino, or m v < 0.23 eV for three degenerate neutrinos. Using the value 
m v = 0.23 eV and multiplying Eq. (13) by 3 for three families we can calculate k. This is 
shown in Figure 2 . The general behavior of k as a function of M* is similar to Figure 1. 
The only point to be noted here is that, because of heavier neutrino mass, the allowed range 
of M* is much higher than the previous case. If we now impose the condition that A < 10 to 
satisfy the perturbative condition, M* becomes too heavy ~ 200 TeV. For such a large value 
of M*, k drops below 10~ 4 for M H = 100 GeV and 5 = 3,, and is essentially unmeasurable. 
The situation is worse for higher values of 5. 

In summary, we see from Figures 1(a), (6) and Figure 2 that k of the order 0(1) can arise 
only when M* ~ 20 — 30 TeV. However, as noted, this implies A > 70, making the theory non 
perturbative. There are two ways by which one can evade this problem. Either one should 



take A ~ 0(1), in that case, M* will be > 100 TeV for m u ~ y Am 2 atm and M* > 200 TeV 
for m v = 0.23 eV. This is perfectly fine, but for such a large value of M#, the invisible decay 
width of Higgs will be negligibly small compared to H — > bb. The second choice suggested 
by [3, 16] is to consider that the singlet neutrino propagates in a sub-space (5„) of the full 



extra-dimension (5) where gravity propagates. Assuming that all extra dimensions are of 
the same size R, in this case, the Dirac mass for the standard model neutrino now becomes 

Thus for 5 U = 5 and 5 = 6 and with M* ~ TeV, A ~ 0(1), we obtain m 2 ~ Am^ m . It has 
been shown by Agashe and Wu [16] that for the above choices of 5 V and 5, the constraints 
on M* from BR{ji — > ej) and n — * eP, fiu decays can be significantly weakened as compared 
to the minimal model, allowing the scale M* ~ TeV. In this case, the maximum value of the 
physical neutrino mass for a given 5 U , S and M* is given by [16]: 



Following Eq. (18), one can invert Eq. (11) to get the neutrino mass parameter m which 
enters in the Higgs decay width. In sub-space, by replacing 5 — > S u and (Mpi/M*) 2 — > 
(Mpi/M*) 2 ^"/ 8 ) one can get the ratio re in terms of the neutrino mass parameter m,6 u ,S, 
and M*. 

« - ^2 hrr hrr ( 19 ) 



3m 2 . \Mj U 4 

We now compute the ratio k as a function of M* keeping 5„ = 5 and 5 = 6. In Figure 3 
(a) and (6) , we show re as a function of M* for = 1.5 x 10~ 3 eV 2 and 4.0 x 10~ 3 eV 2 
respectively. The choice of the Higgs masses is the same as before. Comparing Figure (3) 
with the earlier Figures, it is clear that M* can now be as low as 1 TeV. This value also is 
allowed by other experimental constraints [16]. From Figure 3 (a), one can see that k ~ 1 
only for M H = 200 GeV and at M* = 1 TeV. As M* increases k decreases significantly. 
While in Figure 3 (b), which corresponds to the higher value of the physical neutrino mass, 
m 2 u = 4.0 x 10" 3 eV 2 , re can be ~ 1 even for M H = 150 GeV. But for M H = 100 GeV it 
is always less than 1. For M H = 200 GeV, re can be as large as ~ 10, though for very 
small window in ~ 1 — 1.3 TeV. Similar variation of re with M* is shown in Figure 4 
in the degenerate neutrino mass scenario. In this case, we find that re can be larger than 1 
for M H = 200 GeV. However, as mentioned before, this value of Higgs mass will not serve 
the purpose of looking for invisible decay modes of the Higgs boson. Hence, for practical 
purpose, we should look at values of re for M H up to 150 GeV. It turns out that for M H = 150 
GeV, re ~ 0.8 for M* ~ 2.5 TeV, whereas, for M H = 100 GeV, re can be at most 0.1. 

In these two analysis, we have shown that only in a very small region of parameter space 
can the invisible decay of Higgs be as large as H — > bb decay mode. The main restriction 
comes from the perturbative constraint on the Yukawa coupling A. 



For completeness we also discuss the case of asymmetric dimensions [5, 8]. In this 
scenario neutrinos propagate in a sub-dimensional space of dimension 5 U of size R, whereas 
gravity propagates in space of dimension 5. The extra (5 — 5 V ) has a size r with (r << R). 
In such a scenario Eq. (1) becomes 

M 2 Pl ~ R s ^ M^ +2 r {s - Sv) (20) 

In this case, the Dirac mass for the standard model neutrino becomes 

Af 

m « ^=^= (21) 

To satisfy the constraints from supernova 1987a, the scale 1/R > 10 KeV [7]. For such a 
value of 1/R the mixing angle (~ mR/n) of any KK state with the standard model neutrinos 
becomes negligibly small and m ~ m v . 

One can obtain m 2 u ~ Ara^ tm , for 5„ = 3, A ~ 0(1), M* « 4TeV ( for iT 1 « 10 KeV) 
and M* 8 TeV ( for R^ 1 ps 25 TeV). For <5j, > 3, the neutrino mass wij, is highly suppressed 
as seen from above expression for m (Eq. 21). 

In this scenario, the ratio k as defined in Eq.(13) turned out to be: 



6m 2 b J \M, 



(23) 



To study the invisible Higgs decay in this scenario, we fix A = 1 and take the same values 
of Higgs mass as before. We find k = 0.17,0.6 and 1.3 for M H = 100, 150 and 200 GeV 
respectively for M* = 1.5 TeV. From the Eqs.(22) and (23) it is clear that k decreases as 
5 U , and /or M* increases, while it increases as the higher dimensional Yukawa coupling A 
increases. 

Now we will discuss the observability of this kind of invisible decay of some pseudoscalar 
mesons. 

Neutral Pion decays 

We begin our analysis with the neutral pion. The effective Lagrangian for the process 
7T° — > v L v l R is in the standard model through Z exchange is : 

%G 

£ eff = -^-F n mP-f 5 iy^ (24) 
where, Gf and F n are the Fermi coupling constant and pion decay constant respectively. 



Using this effective Lagrangian we one determine the decay width of pion into vlv 1 r in 
the minimal model: 

We have computed the branching ratio BR(ir° — > v L v l R ) for m 2 = 4.0 x 10~ 3 eV 2 and found 
the result is ~ 1(T 25 for 5 = 3, M* = 200 TeV and A ~ 11. This predicted branching ratio 
is much smaller than the experimental upper limit 8.3 x 10~ 7 at 90% C.L.[17]. 

The branching ratio can be higher by 4 — 5 order of magnitude for lower values of M* ~ 45 
TeV, but such a low value of M* correspond to unacceptably large A ~ 200. In the case of 
degenerate neutrinos, the above branching ratio does not change significantly, so we do not 
present any numerical results here. 

Next, we compute the above decay widths in the sub-space scenario. In this scenario as 
shown earlier one can have A of order one and also M* of about few TeV. With the following 
replacement 5 — > 5 U and (Mpi/M*) 2 — > (Mp^/M*) 2 ^/ 5 ) one can rewrite the decay width 
T (7T° — > vpv l R ) m sub-space model. In this case, we find that for m 2 u = 4.0 x 10 -3 eV 2 , the 
^(tt -> i/ £ pjj) is of the order 10~ 19 for M* ~ 1.4 TeV corresponding to A of 0(1). The 
situation remain unchanged even with the assumption of degenerate neutrino masses. In the 
scenario, where the right-handed neutrinos propagate in the extra- dimensions with largest 
size, the BR(n° — > ul^r) is too small to be observed. 

B meson decays 

In the standard model, the process B — > vv receives contributions from Z-penguin and 
box diagrams, where the dominant contribution comes from intermediate top quark loop. 
Off-shell Z and W exchanges have contributions from would be Goldstone modes that couple 
to right handed neutrinos. The effective Lagrangian for B — > vv for each neutrino is given 
by 

£eff = h^BVV (26) 

where, f\ = 27r S in 2 e w ^t^^td^tb^B^i Gp is the Fermi coupling constants, a is the fine 
structure constant ( at the Z mass scale), Fp is the S-meson decay constant, 9w is the weak 
angle and V t * d V t b are products of CKM matrix elements. The Wilson coefficient Cf^ 1 at the 
leading order is given by: 



where x t = 



Xt + 2 3(x t -2) 1 



(27) 



In models of large extra-dimension, B — > uu gets additional contribution from Higgs 
exchange contribution. The effective flavor changing vertex (bdH°)c&n be obtained from [18] 



G 3 J 2 3 



m b d(l + 75)6 + m d d(l + 75)6 



H° + h.c. (28) 



where, Vij are the elements of the Kobayashi-Maskawa matrix and mi are the corresponding 
quark masses flowing in the loop. 

The effective Lagrangian for the process B — > uu retaining only the top quark contribu- 
tion is: 

C eff = f 2 *Bi&L (29) 



where, f 2 = %\{^G F f' 2t -^^V td V tb ^^r\ M H is the standard model Higgs mass and 
v is the standard model vacuum expectation value. The decay width B — > uip l R (after 
summing over all KK modes of right handed neutrinos below ms) can be obtained after 
adding contributions from Z-penguin and box diagrams and the Higgs mediated diagrams 
together 

TV D j-.i \ _ f total ( m B\ 5 ( M Pl\ 2 m /o nl 

F(B - u L u R ) - — ^— ) ^_ ) m B (30) 

where, f total — fi + /2- To compute the branching ratio we have taken following input 
parameters m b = 4.2 GeV, m B = 5.279 GeV, m t = 175 GeV, v = 246 GeV, V td = 0.006, 
V tb = 1.0, F B = 0.180 GeV. We then obtain the BR(B — > u L u R ) for case of degenerate 
neutrinos with m u = 0.23 eV ( which also corresponds to the heaviest neutrino) and find 
that the branching ratio varies between ~ 10~ 12 — 10~ 14 for M* ~ 50 — 120 TeV, 5 = 3 and 
A ~ 80 — 20. This branching ratio is too small to be observed at any present B factories. 
Even in the sub-space scenario, the branching ratio does not get any significant enhancement, 
irrespective of the different neutrino masses, assumptions considered previously. In the 
scenario, where the right-handed neutrinos propagate in the extra- dimensions with largest 
size, the BR(B — > ulu r ) is too small to be observed. 

Conclusions 

In this analysis we have studied the possible enhancement of invisible decay widths of the 
standard model Higgs boson and other pseudoscalar mesons in the model of singlet neutrinos 
in extra-dimensions. In the case of Higgs boson decay we have found that in certain range 
of extra-dimension parameter space, the branching ratio of Higgs into invisible mode can be 
> BR(H — > bb). Unfortunately, the higher dimensional Yukawa coupling A takes on large 



(> 100) values in that parameter space. For A < 10, the H — > vv rate is a tiny fraction of 
the H — > bb rate. 

We have also studied the invisible decay rate of Higgs in the scenario, where right-handed 
neutrinos are in sub-space (5 V < 5), which is the modification of the minimal model, required 
to keep A < 0(1). It has been shown that to have a consistent model allowed by different 
experimental constraint , one should have 5 U = 5 and 5 = 6. In this scenario, the invisible 
decay rate of Higgs can compete to that of H — > bb, though for a very small range of 
~ 1 — 1.3 TeV. In the scenario, where the right-handed neutrinos propagate in the 
extra-dimensions with largest size, the invisible decay of Higgs can be as large as H — > bb 
for Mu = 200 GeV, 5 U = 3, A = 1 and M* = 1.5 TeV. However, for lower values of Mn 
(= 100 and 150 GeV) the invisible decay rate is smaller than H — > bb for the above set of 
parameters. 

We have also studied the decay rates 7r° — > vl^r and B — > vl^r in all these scenarios. 
Unfortunately in both of these decays the new effects are negligibly small to be measured. 
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Figure 1: Variation of k (defined in the text) with M* for three values of extra dimensions 
5 = 3,4 and 5. The solid, broken and dotted lines correspond to Mh = 100, 150 and 200 
GeV respectively. Figures (a) and (b) correspond to = 1.5 x lCT 3 eV 2 and 4.0 x 10~ 3 eV 2 
respectively. 




Figure 2: Variation of k with M* for three values of extra dimensions 5 = 3,4 and 5. The 
solid, broken and dotted lines correspond to M H = 100, 150 and 200 GeV respectively. We 
have fixed m u = 0.23 eV. 




Figure 3: Variation of k with M* for number of sub-space extra dimensions S v = 5, number 
of extra-dimensions 5 = 6 and M H = 100, 150 and 200 GeV. Figures (a) and (b) correspond 
to ml = 1.5 x 10~ 3 eV 2 and 4.0 x 10~ 3 eV 2 respectively. 
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Figure 4: Variation of k as a function of M* in sub-space extra-dimensions 5 U = 5, number 
of extra-dimensions 5 = 6 and M# = 100, 150 and 200 GeV. We have fixed m v = 0.23 eV. 



